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Discrete Data Assimilation in the Lorenz 
and 2D Navier— Stokes Equations 

Kevin Hayden 1 , Eric Olson 2 and Edriss S. Titi 3,4,5 

Consider a continuous dynamical system for which partial information 
about its current state is observed at a sequence of discrete times. Dis- 
crete data assimilation inserts these observational measurements of the 
reference dynamical system into an approximate solution by means of 
an impulsive forcing. In this way the approximating solution is coupled 
to the reference solution at a discrete sequence of points in time. This 
paper studies discrete data assimilation for the Lorenz equations and the 
incompressible two-dimensional Navier-Stokes equations. In both cases 
we obtain bounds on the time interval h between subsequent observa- 
tions which guarantee the convergence of the approximating solution 
obtained by discrete data assimilation to the reference solution. 

1. Introduction 

In [11] and [12] Olson and Titi studied the number of determining modes for continuous 
data assimilation for the incompressible two-dimensional Navier-Stokes equations. As in 
those papers, the motivating problem for our work is the initialization of weather forecast- 
ing models using near continuous in time measurement data obtained, for example, from 
satellite imaging. In this work, rather than making the idealization that the measurement 
data is continuous in time, we focus on the case where the measurement data is taken at 
a sequence of discrete times t n . 

If t n = hn and h is small then discrete data assimilation can be viewed as near 
continuous. One expects that if the approximating solution obtained by continuous data 
assimilation converges to the reference solution then the approximating solution obtained 
by discrete data assimilation for small h will also converge to the reference solution. Note, 
however, that near continuous observational data is mathematically quite different from 
continuous data. If the observations are known continuously in time on some interval then 
mathematically the n-th time derivatives may be calculated on that same interval for all 
values of n. It is possible that this derivative information could lead to a reconstruction 
of the reference solution in cases where near continuous measurement information might 
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not. For example, Wingard [18] shows for the Lorenz equations that knowing X and its 
time derivatives at a single point in time can be used to recover both Y and Z. 

In this paper we present a technique to prove a discrete in time determining mode 
result in the specific context of creating an approximating solution that asymptotically 
converges to a reference solution. In addition to providing a more realistic framework in 
which to study near continuous data assimilation, our work can be seen as a discrete in 
time extension of the theory of determining modes developed by Foias and Prodi [5] and 
further refined by Jones and Titi in [7] and [8]. Unless otherwise noted, we shall assume 
the dynamics governing the evolution of the reference solution admit a global attractor 
and that the reference solution lies on that global attractor. This assumption is made 
for simplicity, as our analysis actually depends only on the existence of an absorbing ball 
which contains the reference solution forward in time. 

We begin our discussion with the easier case of the Lorenz equations to provide insight 
and illustrate the methods we will use for the incompressible two-dimensional Navier- 
Stokes equations. General studies for the synchronization of discrete in time coupled 
systems for the Lorenz equations were produced by Yang, Yang and Yang [20] and Wu, 
Lu, Wang and Liu [19]. Note, however, that the matrix corresponding to the observational 
measurements studied here does not have a suitable spectral radius to apply the theorems 
of their work. 

The method of discrete data assimilation can be described mathematically as follows. 
Let U be a solution lying on the global attractor of a dissipative continuous dynamical 
system with initial condition Uq at time to- Let S be the continuous semigroup defined by 
U(t) — S(t, to, Uq). Represent the observational measurements of the reference solution U 
at time t n by PU (t n ), where P is a finite-rank orthogonal projection and t n is an increasing 
sequence in time. Discrete data assimilation inserts the observational measurements into 
an approximate solution u as the approximate solution is integrated in time. In particular, 
let uq = r) + PU(to) and u n+ i = QS(t n+ i,t n ,u n ) + PU(t n+ i) for n = 0,1,..., where 
Pi] = and Q = I — P. Here r\ corresponds to an initial guess for the part of the reference 
solution QU (to) that can not be measured. The approximating solution u obtained by 
discrete data assimilation is defined to be the piecewise continuous in time function 



Our goal is to find conditions on P, t n and rj which guarantee that the approximating 
solution u converges to the reference solution U as t — > oo. 

For the Lorenz system the reference solution is a three dimensional vector consisting 
of the components X, Y and Z whose evolution is governed by the coupled system of three 
ordinary differential equations 



where a = 10, b = 8/3 and r = 28. We shall assume that the reference solution lies on the 
global attractor. 



u{t) = S(t,t n ,u n ) for te[t n ,t n+1 ). 



(1.1) 




bZ + XY- b{r + a) 



(1.2) 
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We take the observational measurements of the reference solution to be the values of 
the variable X at the times t n . These observations of X are used to create an approximating 
solution whose components are x, y and z. Note that x(t n ) = X(t n ) where y(t) and z(t) 
are continuous at t = t n and x, y and z satisfy 

x = —ax + ay 

y = —ax — y — xz (1-3) 
i = — bz + xy — b(r + a) 

on each interval [t n , t n _i) for n = 0, 1, 2, ... . 

For simplicity we take t n = hn for some fixed h > 0. Numerical experiments for the 
Lorenz system done by Hayden [6] indicate that the approximating solution converges to 
the reference solution as t — )■ oo for values of h as large as 0.175. Figure 1.1 shows the 
convergence of the approximating solution to the reference solution when h = 0.1. 



4 







-4 



-6 



i \ 
V 

t k 

A \ 
\ \ 


1 1 1 

X x 
Y-y 

- Z-z 




l \ \ . ' s V / 

I \ V: ! i \ A * 

\ \; - 

V si ' 
\ ft 

V A ' 

- 

1 1 1 1 1 1 1 1 1 



;| 1 , , , 1 , 1 , , 1 

0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 



Figure 1.1. Convergence of the approximating solution to the reference 
solution for the Lorenz system when h = 0.1. 

Note again that the variables y and z in the approximating solution are continuous in time, 
in particular at t = t n , whereas x is discontinuous. The discontinuities in x are the result 
of the assimilation of the observations of X at each time t n . 

Our main result for the Lorenz equations is an analytic proof of Theorem 2.5 which 
shows there exists t* > depending on a, /3 and r such that for any h G (0,t*] the 
approximating solution obtained by discrete data assimilation of measurements of the X 
variable at times t n = hn, as described in (1.3) above, converges to the reference solution 
as t — > oo. 
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The second part of this paper focuses on the incompressible two-dimensional Navier- 
Stokes equations with L-periodic boundary conditions 

f -vAU + (U-V)U + VP = f (14) 

V -U = 

where v is the kinematic viscosity and / is a time independent body forcing. In this case 
U can be expressed in terms of the Fourier series 

U=^2ll k e ik ' x where J = j -j-( n i, n 2) ■ e Z and {n u n 2 ) ^ (0, 0) 1 
kej ^ ' 

where Uk • k = for all k e J" and £/& = £7* fc . Again, we assume that the reference solution 
lies on the global attractor. 

Let P\ be the orthogonal projection defined by 

P X U = U ^' X - (1-5) 

k 2 <\ 

We take the observational measurements of the reference solution to be the values of P\U 
at the times t n . Note that A -1 / 2 represents the smallest length scale of the fluid which can 
be observed — the resolution of the presumed measuring equipment or distance between 
measuring stations. 

In the context of continuous data assimilation and determining modes, the rank of 
the smallest projection P\ such that the approximating solution converges to the reference 
solution is called the number of determining modes. For discrete and near continuous data 
assimilation, the parameter A also depends on the interval of time between the observational 
measurements. Our main results for the incompressible two-dimensional Navier-Stokes 
equations are Corollary 3.10 and Corollary 3.11. 

Corollary 3.10 shows for any time interval h > there exists A large enough so that 
the approximating solution obtained by discrete data assimilation of the measurements 
P\U(t n ) where t n = hn will converge to the reference solution. This means that increasing 
the resolution of the measurements can compensate for a large time interval between 
subsequent observations. Corollary 3.11 shows there is a dimensionless constant C such 
that if 

then there exists h > depending only on |/|, v, O and A small enough such that the 
approximating solution obtained by discrete data assimilation with initial guess rj = 
converges to the reference solution. This means that decreasing the time interval between 
subsequent observations can compensate for a low resolution provided the resolution meets 
a minimum standard. 
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2. Lorenz Equations 

Following Foias, Jolly, Kukavica and Titi [4] we write the Lorenz system (1.2) as 



dU 
~dt 



+ AU + B(U,U) = f 



(2.1) 



where 
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and again a = 10, b = 8/3 and r = 28. One reason for writing the Lorenz equations 
in this way is to make the similarities and differences in the proofs from this section and 
following section on the Navier-Stokes equations more transparent. Despite the notational 
similarities there should be no trouble distinguishing the results on the Lorenz equation 
that apply only to the Lorenz equations from the results on the Navier-Stokes equations 
which apply only to the Navier-Stokes equations. We start with some definitions and facts 
that are easy to verify. 

Definition 2.1. \U\ = y/(U, U) = VX 2 + Y 2 + Z 2 . 

The facts below may be deduced from the preceding definitions and are listed here for 
reference. First, we state an estimate involving the linear term AU that plays the same 
role in our treatment of the Lorenz system that Theorem 3.2 plays in our later treatment 
of the incompressible two-dimensional Navier-Stokes equations: 



(AU, U) = aX 2 + Y 2 + bZ 2 > \U\ 



(2.2) 



Next we state some algebraic identities analogous to the orthogonality relations (3.3) and 
(3.4) for the nonlinear term in the two-dimensional Navier-Stokes equations: 



(B(U,U),U) = and B(U,U) = B(U,U). 
We will also use an estimate on B which is similar to Theorem 3.5: 

\B(U,U)\ < 2~ 1 \U\\U\. 



(2.3) 



(2.4) 



The following bound on the global attractor was shown in Temam [16] on page 33, 
see also Foias, Constantin and Temam [3]. 

Theorem 2.2. Let U be a trajectory that lies on the global attractor of (2.1). Then 
\U(t)\ 2 < K for all £ G R where 

b 2 (r + a) 2 

K= 4(6-1) • (2 - 5) 
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Before proceeding, we state Young's inequality which will be used throughout the 
remainder of this work. 

Theorem 2.3. Let 1/p + l/q=l then \xy\ < \x\ p /p + \y\ q /q. 

We are now ready to begin our study of discrete data assimilation for the Lorenz 
equations. Our measurements will consist of the X variable at the times t n . Therefore, we 
define the orthogonal projection P as 



1 





and Q = I - P 




1 
1 



Let U be a solution of the Lorenz equations lying on the global attractor and u the 
approximating solution given by (1.1) where S is the semigroup generated by (2.1). Since 
U and u both satisfy (2.1) on the time interval [t n ,t n +i) setting 5 = U — u we obtain 

— + uA5 + B(U,U)-B(u,u) = for te[t n ,t n+1 ) 
at 

or after some algebra that 

— + AS + B(U,5) + B(5,U)-B(S,5) = for te[t n ,t n+1 ). (2.6) 



Lemma 2.4. There exists (3 > given by (2.8) below, depending on a, b and r, such that 

\5(t)\ 2 <\5(t n )\ 2 e^-^Uorte[t n ,t n+1 ). 

Proof: Take the inner product of S with (2.6) and apply (2.2) and (2.3) to obtain 



ld\S\ : 



2 dt 



+ \5\ 2 + 2(B(U,5),5) < 0. 



Estimating using (2.4) and Theorem 2.2 gives 

2\(B(U,S),S)\ < \U\\5\ 2 < K l ' 2 \5\ 2 . 

Therefore 



d\5[ 
dt 



where 



P = 2(K 1 ' 2 - 1) 
and K is defined in (2.5). Integrating from t n to t yields 

\S(t)\ 2 < |<5(t n )| 2 e^-H 

which finishes the proof. 



(2.7) 



(2.8) 



□ 
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Theorem 2.5. Let U be a solution of the the Lorenz equations (2.1) lying on the global 
attractor. Then, there exists t* > depending only on a, b and r such that for any 
h G (0,t*] the approximating solution u given by (1.1), see also (1.3), with t n = hn 
converges to U as t — > oo. 

Proof: Take the inner product of PS with (2.6) to obtain 



1 d\PS\' 

2 dt 



+ (AS, PS) = 0. 



Let 





' -a 0" 




'5y~ 
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Then \P5\ 2 = 8 2 x and 



(AS, PS) = (PAS, PS) = (APS, PS) + (w, PS) = aS 2 x - aS x Sy. 
Therefore (2.9) becomes 



1 d8 x 9 a r9 a r9 _ a r9 a . „. 



where we have applied Theorem 2.3 with p = q = 2. Using Lemma 2.4 we obtain 

d\PS\ 2 



dt 



+ a\PS\ 2 < aS 2 < a\S(t n )\ 2 e 



•2 p P(t-t n ) 



Multiplying by the integrating factor e a ^ and integrating from t n to t gives 



\PS(t)\ 2 < 



<r\5(tn)\ 

P + <7 



■0 



,/3(t-t n ) _ -a(t-t n ) 



) 



where we have used the fact that PS(t n ) = 0. 

A finer analysis of the nonlinear term appearing in (2.7) gives 



Therefore 



2(B(U,S),S) = S X (ZS Y -YS Z ). 



\2(B(U,S),S)\ < \U\\PS\\S\ < K 1/2 |P5||5| < \k\PS\ 2 + ^-\S\ 2 . 



Substituting (2.10) and (2.11) into (2.7) yields 



d\S\'< 
~dt 



+ \S\ 2 < K 



°\S(t n )\< 
P + a 



-<j(t-t n ) 



(2.9) 



(2.10) 



(2.11) 



(2.12) 
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Multiply (2.12) by *") and integrate from t n to t to obtain 

m)\ 2 <M(t-t n )\d(t n )\ 2 

where 

M(r) = e~ T (l + ^ (e^ +1 ) s - e^-^ds) (2.13) 

is a function that doesn't depend on t n . Note that M(0) = 1. Differentiating yields 

M'(r) = -M(r) + 4^—{J T ~ e ~ aT )- 
p + a 

Therefore M'(0) = -1. It follows that there is t* > such that 

M(h) = 1 + / M'(s)ds < 1 

for all h G (0,**]. 

Next fix h G (0,**] and let 7 = M(/i) < 1. Then 

|5(t n+1 )| 2 = |Q5(t n+1 )| 2 = lim |Q5(t)| 2 < lim \5(t)\ 2 

< lim M(t-t n )\5(t n )\ 2 = M(h)\5(t n )\ 2 = 1 \5(t n )\ 2 

implies by induction that 

Wn)\ 2 < i n Wo)\ 2 = i n \QU(t ) - v \ 2 < i n R 

where R = 2(K + \i]\ 2 ). Now let t > and choose n so that t G [t n , t n+ \). Then n — >■ 00 as 
t — )> 00 and therefore 

|5(t)| 2 < M(t - t n )|5(t n )| 2 < |5(t n )| 2 < 7 n # (2.14) 

shows that the approximating solution converges to the reference solution as t — > 00. □ 
Corollary 2.6. If a = 10, b = 8/3 and r = 28, then £* « 0.000129. 

Proof: The result follows using the value of iiT from Theorem 2.2 and choosing t* slightly 
less than the value of t such that M(t) = 1. □ 

Before proving the final result in this section, we extend Theorem 2.5 by proving that 
the interval of time between t n +i and t n need not be exactly h for every n G N. 

Corollary 2.7. Let t* be the bound given in Theorem 2.5. Suppose £ n +i —t n <t* where 
t n — >■ 00 as n — > 00. Then the approximating solution u given by (1.1) converges to the 
reference solution U of the Lorenz equations (2.1) as t — > 00. 
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Proof: Let h n = t n — t n -\. If there exists e > such that the set /C = { k : h k > e } is 
infinite then M(h k ) < max { M (s) : s G [e, t*] } < 1 for k G K implies that 

n 

\S(t n )\ 2 <Y[M(h k )R^0 as n ->• oo. (2.15) 

fc=i 



Otherwise, /i n ->■ as n ->■ oo. By Taylor's theorem \M(h n ) - M(0) - ft n M' (0)| < C/i 2 
where C = | max { \M"(s)\ : s G [0,**] }. Choose N so large that Ch k - 1 < -1/2 for 
k > N. Since M(0) = 1, M'(O) = -1 and logx < a; - 1 for x > it follows that 

n n n 

£ logM(fe fc ) < £ (M(^ fc ) - 1) < £ ( M (°) + ^fc^(O) + Ch 2 k ~ 1) 

fc=JV fc=JV k=N 

n 1 n 

hk(Chk — 1) < — £^ /ifc — > — oo as n — > oo. 



Thus 



2 

k=N k=N 



n n 

j J M(h k )R = exp ( log i? + £ log M(h k f) ^0 as 



n — > oo. 

fc=i fc=i 
The proof now finishes as in Theorem 2.5. □ 

Our final result on the Lorenz equations shows that the approximating solution is 
bounded for any updating time interval of h. If h < t* then the approximating solution 
converges to the reference solution, and since the reference solution is bounded, then 
the approximating solution will also be bounded. In the case where the approximating 
solution does not converge to the reference solution then the following result shows that 
the approximating solution is still bounded. 

Theorem 2.8. Let U be a trajectory that lies on the global attractor of (2.1). The 
approximating solution u given by (1.1) where t n = nh with h > is bounded. Namely, 
there is a constant M x that depends only on rj, a, b and r such that \u(t) | 2 < Mi/(1 — e~ h ) 
for all £ > and h > 0. 

Proof: Taking inner product of (2.1) with u and using (2.2) followed by Young's inequality 
we obtain 

^ + l«P<(/,»)<>^M I 

and consequently 

/ '" l2 +l«| 2 <l/| 2 (2-16) 



d\u\ 2 



dt 

for t G [t n ,t n +i). Gronwall's inequality then implies 

»(t)| 2 < |«„| 2 e-<*-*-> + |/| 2 (1 - e-<*-*-)) for te[t n ,t n+1 ). 



Defining 



u =r] and u n+ i = lim u(t) = S(t n+ i, t n , u n ) 
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so that u n — Qu n + PU(t n ) we obtain 

K +1 | 2 < lim (| Un | 2 e -(^) + |/| 2 (l- e -(^))) 
tyt n+1 

<\Un\ 2 l+\f\ 2 (l-l) 

where 7 = e~ h . Since 

\u n ? = \Qu n \ 2 + \PU(t n )\ 2 < \u n \ 2 + \U(t n )\ 2 < \u n \ 2 + K 

we obtain 

\u n+1 \ 2 <\u n \ 2 1 + C 1 (2.17) 

where d = K 7 + \f\ 2 (l - 7). 

Induction on (2.17) and summing the series yields 

\u n \ 2 < |uo|V + Ci(l + 7 + • • • + 7 n_1 ) = H 2 l n + Ci—^. 

1-7 

Given t > choose n so that t e [t n , t n +i)- Then t = t n + a where a G [0, /i) and 
\u(t)\ 2 < \u n \ 2 e~ a + \f\ 2 (l - e~ a ) < \u n \ 2 + |/| 2 

<\u n \ 2 + K+ |/| 2 < |r,| V + Ci-^- +K + |/| 2 < -— i- 

1 — 7 1 — 7 

where Mi = |?7| 2 + Ci + K+ |/| 2 . □ 

If we take rj = as our initial guess for QU(to) when forming the approximating 
solution u then the constant Mi depends only on a, b and r. In either case we obtain the 
asymptotic bound 

lim sup |w(t)| 2 < 



,2 . ^2 



£—►00 



1 -e- fc ' 



where the constant M 2 depends only on a, 6 and r. 

Theorem 2.8 can be improved using the exact structure of (Au,u) and (f,u). In 
particular, we obtain using Theorem 2.3 that 

(Au, u) - (/, it) = ax 2 + y 2 + bz 2 + zb(r + a) 

= ax 2 + y 2 + z 2 + (b- l)z 2 + zb(r + a) 

> 2^ 2^ 2 b 2 {r + a) 2 , 2 |/| 2 

> ax +y +z — — — > \u\ 



4(6-1) - 1 1 4(6-1) 

which improves the bound to |w(t)| 2 < M%/{1 — e~ 2h ). Although an improvement, this 
bound on u still tends to infinity as h tends to zero. As mentioned earlier, a bound uniform 
in h can be obtained by combining Theorem 2.5 with Theorem 2.8. 

Corollary 2.9. There exists a bound M4 depending only on 77, a, 6 and r such that the 
approximate solution u obtained with t n = hn is bounded by M4 for any h > 0. 
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Proof: Let t* be given as in Theorem 2.5. If h < t* then (2.14) implies 
\u(t)\ = \U(t) - 5{t)\ < \U(t)\ + \5{t)\ < K 1 ! 2 + R 1 / 2 . 
If h > t* then Theorem 2.8 implies 

|«(f)| < Ml /2 /(l - e- h ) 1 ' 2 < Ml /2 /(l - e" r 

Taking 

M 4 = max^ 1 / 2 + R 1 / 2 , Ml /2 /{l - e" r ) 1 / 2 } 

yields a constant that depends only on rj, a, b and r. □ 

As mentioned before, if we take 1] = or consider asymptotic bounds on u as t — > oo 
then the dependency on r\ can be removed from these bounds. 



3. Navier— Stokes Equations 

This section contains results for the two-dimensional incompressible Navier-Stokes equa- 
tions that are similar to the results proved in the previous section for the Lorenz equations. 
Let O = [0, L] x [0, L], V be the space of divergence- free vector- valued L-periodic trigono- 
metric polynomials from O into R 2 with J n u = 0, H be the closure of V with respect 
to the L 2 norm and Ph be the L 2 -orthogonal projection Ph'-L 2 (Q) — > H, referred to 
as the Leray-Helmholtz projector. Following the notations of Constantin and Foias [2], 
Robinson [15] and Temam [17] we write the incompressible two-dimensional Navier-Stokes 
equations (1.4) as 

— + isAu + B(u,u) = f (3.1) 

LLC 

where 

Au = -P H Au and B(u, v) = P H [(u ■ V)u] . (3.2) 

Note that we have assumed / G H so that Pnf = f ■ Also notice that A = —A in the 
periodic case. 

Let V be the closure of V with respect to the H 1 norm and V(A) = H n H 2 (Q) be 
the domain of A. We may define norms on H, V and V(A) which are equivalent to the 
L 2 , H 1 and H 2 norms respectively by 



\u\ 



L 2 |w/c| 2 , \\u\\ = L 2 k 2 \uk\ 2 and \Au\ = L 2 fc 4 | 
kej kej kej 



Here u has been expressed in terms of the Fourier series 



u = ^ Uke lk ' x where J = j -^-(ni, 77-2) : ni e Z and (ni, 712) 7^ 1 . 
kej ^ ' 



The mathematical theory proving the existence and uniqueness of strong solutions to 
the two-dimensional incompressible Navier-Stokes equations (3.1) with initial data in V 
may be found for example in [2], [15] or [17]. Specifically we have 
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Theorem 3.1. Let uq G V and / G H. Then (3.1) has unique strong solutions that satisfy 

ua»((0,T);F)ni 2 ((0,T);^)) and ^ G L 2 ((0, T); tf) 

for any T > 0. Furthermore, this solution is in C([0,T]; V') and depends continuously on 
the initial data uq in the V norm. 

Let S{t,to,uo) for t > t denote the unique strong solution of (3.1) given by 
Theorem 3.1 with initial condition uq at time to- Let Ai = (2n) 2 /L 2 be the smallest 
eigenvalue of A on T>(A) and define P\ as in (1.5). Further define Q\ = I — P\ to be the 
orthogonal complement of P\. The Poincare inequalities for u, P\u and Q\u on O can 
now be summarized as 

Theorem 3.2. Given P\, Q\ and Ai as defined above then 

\u\ < X~ 1/2 \\u\\, \Q\u\ < \- 1/2 \\Q x u\\, \\P\u\\ < \ 1/2 \P x u\ 

and 

IMI < \i 1/2 \Au\, \\Q\u\\ < X~ 1/2 \AQ x u\, \AP x u\\ < A 1/2 ||P A w||. 
provided the norms exist and are finite. 
Theorem 3.3. If u G V(A) then 

\Au\ l ' A < 2 1 /8(A 1 /8|| W ||V4 + |Q A Au| 1 / 4 ). 

Proof: Since P\ and Q\ are orthogonal then 

\Au\ 2 = \P x Au\ 2 + \Q x Au\ 2 < X\\Pxu\\ 2 + \Q x Au\ 2 < 2 max{A||P A w|| 2 , |Q A ^| 2 }. 
Therefore 

{An] 1 '* < 2 1 /8 ma x{A 1 /8||P AW ||V4 7 |Q A ^ u |l/4} < 2 1 /8(A 1 /8|| U ||V4 + IQ^) 1 / 4 ), 

which completes the proof. □ 

Let us now recall some algebraic properties of the non-linear term that can also be 
found in [2], [15] or [17] that play an important role in our analysis. They are 

(B(u,v),w)=(B(u,w),v), (3.3) 

(B(u,v),v) = and (B(u, u), Au) = 0. (3.4) 
We now move on to some inequalities which we shall refer to later. 
Theorem 3.4. There exists a dimensionless constant c such that if u G V then 

IMlls/3 < C|w| 3/2 |M| 1/2 , \\u\\ 2 L 4. < C\u\\\u\\, \\u\\ 2 L s < C|w| 1/2 ||w|| 3/2 
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and if u E V(A) then 

NU~ < c a;; 1/8 hI 3/4 W 1/4 - 

The first three inequalities above may be obtained from the Sobolev inequalities fol- 
lowed by interpolation, see for example (6.2) and (6.7) in [2]. The inequality bounding L 4 
is sometimes referred to as Ladyzhenskaya's inequality and appears as Lemma 1 on page 8 
of Ladyzhenskaya [10]. The last inequality is a form of Agmon's inequality which appears 
as (2.23) on page 11 in Temam [17]. These inequalities may be used along with Holder's 
inequality to estimate the nonlinear term. 

Theorem 3.5. Let u,v,w <EV then 

\(B(u, v),w)\ < c|-u| 1/2 |M| 1/2 || v || \w 

l 1/2 IMI 1/2 

if further v eV(A) then 

\(B(u,v),w)\ < cA-^iMHMi 3 / 4 ^! 1 / 4 ]™! 

\(B(v,u),w)\ < c\~ 1/8 \\u\\\\vf /4 \Av\ 1/4 \w\ 

where c is the constant appearing in Theorem 3.4. 

We consider a reference solution U to the incompressible two-dimensional Navier- 
Stokes equations and the approximating solution u given by (1.1) where S is the semigroup 
generated by (3.1). As with the Lorenz equations, we shall assume that the reference 
solution lies on the global attractor. As proved in Jones and Titi [7] we have 

Theorem 3.6. A solution U that lies on the global attractor of (3.1) satisfies the bound 
\\U\\ 2 < K where 

Ifl 2 

K = ^. (3.5) 

A similar result also appears in [2] for establishing estimates on the global attractor. 
As mentioned in the introduction, any reference solution U which satisfies a bound such 
as (3.5) forward in time is suitable for our analysis, whether that solution is on the attractor 
or not. However, for simplicity we continue to assume U lies on the global attractor. 

Now setting 5 = U — u and following the same algebra as for the Lorenz equations we 
arrive at 

— + uA6 + B(U, 5) + B(6, U) - B{5, 5) = 0, (3.6) 

an equation that looks like (2.6) where A and B have been given the meanings in (3.2). 

Before starting with the proof of our main result we begin with a simpler theorem 
that finds values of A large enough for the approximating solution u to start converging to 
a reference solution U lying on the attractor but does not provide the uniform bound on h 
necessary to maintain convergence as t — > oo. 

Theorem 3.7. If A > c 2 |/| 2 /(Aiz/ 4 ) there exists t > t n such that \6(t)\ < \S(t n )\. 
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Proof: Multiply (3.6) by S and integrate over O. The first of the orthogonality relation- 
ships given in (3.4) yields 



ld\5 



- + u\\S\\ 2 = -(B(8,U),8). (3-7) 



2 dt 

Now applying Theorem 3.5 followed by Young's and Theorem 3.6 we obtain 

^ + H|5|| 2 <c|5|||5||||t/||<^||5|| 2 + ^H 2 

or 

dt ~ v 

Integrating from t n to t where t G [t n , t n+ i) results in 

\5(t)\ 2 - \5(t n )\ 2 < f (^l^)| 2 - V\\6(8)f)d8. 

Therefore 

\S(t)\ 2 <M n (t-t n )\5(t n )\ 2 

where 

Mn(r) = 1 + £ ^ (^f\S(tn + s)\ 2 - u\\5(t n + S )f)ds. 

Note M n (r) is a differentiable function with M n (0) = 1. 
Differentiating yields 

M , ( , c 2 K \5(t n + r)\ 2 J5(t n + r)\\ 2 

M n{ T ) = — \ XU M2 V ISU M2 ' 



\S(tn)\ 2 |<*(*„), 

By Theorem 3.2 and the hypothesis A > c 2 K/u 2 we obtain 

<(0) = ff* _ , M = i* _ , iWf < _ „ A < p. 

i/ |5(t n )| 2 z/ |5(t n )| 2 1/ 

Therefore there exists r > such that 

M 



(r) = 1 + / M'(s)cfe < 1. 



It follows that for t = t n + r that \6(t)\ < \S(t n )\. □ 

It should be pointed out that the value t in the above theorem depends on M n (r) and 
thus on S(s) for s > t n . Unfortunately, this provides us with no way of knowing whether 
the family of functions M n is equicontinuous or not. This means Theorem 3.7 can not 
be used to provide a uniform bound t* on the discrete measurement time interval h that 
ensures the approximate solution will converge to the reference solution as t — > oo. In 
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Theorem 3.9 below we circumvent this issue at the expense of a more stringent condition 
on the minimum size of A. 

Note that the bound on A given in Theorem 3.7 above is the same as the bound given 
by Jones and Titi in [7] on the number of determining modes for the incompressible two- 
dimensional Navier-Stokes equations. This bound was later improved by Jones and Titi 
in [8] using estimates on the time averages of the term ||Att|| 2 . However, in the case of 
discrete assimilation the same technique did not achieve similar improvements. 

We commence with the detailed analysis that allows us to prove our main result for 
the incompressible two-dimensional Navier-Stokes equations. The first theorem we prove 
is an analog of Lemma 2.4. 

Theorem 3.8. There exists /3 > depending on |/|, O and v such that the solution 5 to 
(3.6) satisfies \\d(t)\\ 2 < ||<?(t n )|| for t G [t n ,t n+1 ). 

Proof: Multiply (3.6) by A5 and integrate over O to obtain 

By Theorem 3.5 we have that 

|(B(5,t/),^)|<cAr 1/8 ||5|| 3 / 4 ||C/|||^| 5 / 4 

and 

\(B(U,S),AS)\ < cA~ 1/8 ||5|| 3/4 ||?7|||A5| 5/4 . 
Therefore applying Theorem 2.3 with p = 8/3 and q = 8/5 we obtain 



|2 



l d\\S n 
2 dt + ^ 



A5\ 2 <2c\- 1/8 \\5\\ 3 / 4 \\U\\\A5\ 5 / 4 

< Ciz/- 5 / 3 A7 1/3 ||5|| 2 ||t/|| 8/3 + v\A5\'' 



where C\ is a dimensionless constant related to c. Hence 

d\\bf 



dt 



where ft = 2C\v 5 / 3 A 1 1 / 3 i^ 4 / 3 . Integrating in time from t n to t and applying Gronwall's 
inequality yields \\5(t)\\ 2 < ||<y(t„)|| 2 e^(*-*») for t G [t n ,t n+1 ). □ 

We are now ready to prove the main results of this paper. 

Theorem 3.9. If 

A> 9^ 2c|/|/(AV 2 z,)+c||^ )|| x8/3 ^ 

then there exists t* > depending only on K, \\5(to)\\, u, O and A such that for any 
h G (0, t*] the approximating solution u given by (1.1) with t n = hn converges to the 
reference solution U of (3.1) as t — > oo. 
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Proof: The proof is similar to the proof of Theorem 2.5 for the Lorenz equations with the 
addition that we first use induction to show the bound R = \\S(to)\\ 2 on the difference of 
the initial conditions ensures that ||5(t n )|| 2 < R holds for each each t n . 
Define 

^^MA^'f^^+MA-f 3 ^ (3 - 9) 

where 

C 2 = c 2 (2 5 / 4 ), C 3 = 3c 8 / 3 (5 5 / 3 )(2- 10 / 3 ) (3.10) 

and 

gi (r) = e^{R 1 ' 2 e^/ 2 + 2K 1 ' 2 )\ g 2 (r) = ef* {R 1 ' 2 e^' 2 + 2K 1 / 2 f\ 
Further define 

M(r) = e~ uXT (l + g(s)e uXs d s y (3.11) 
Note that M(0) = 1. Differentiating yields 

M'(r) = -u\M(r)+g(T). 

Therefore 

M'(0) = -v\ + C2(^) 1/ \R 1/2 + 2K 1 ' 2 ) 2 + C 3 [-^-) 1/3 {R 1 / 2 + 2K 1 ' 2 f\ 

Taking C 4 = 9c 8 / 3 > max{(2C 2 ) 4 / 3 , 2C 3 } = 3c 8 / 3 (5 5 / 3 )(2" 7 / 3 ) we find that M'(0) < is 
guaranteed when 

1 ) ■ 

Now choose t* > such that M(h) < 1 for all h G (0,£*]. Note that t* only depends 
on K, \\S(to)\\, u, O and A. We now show the approximating solution u given by (1.1) with 
t n = hn converges to the reference solution U of (3.1) as t — > oo. 

For induction on n we suppose that ||5(t n )|| < R. In the case n = the induction 
hypothesis is true by definition. Take inner product of AQ\5 with (3.6) to obtain 

\j t WQM 2 + "\AQxS\ 2 <B 1 + B 2 + B 3 

where 

B 1 = \(B(S,U),AQ X S)\, B 2 = \(B(U,S),AQ X S)\ and B 3 = \(B(S, S), AQ X S)\. 
Estimate using Theorem 3.5 followed by Theorem 3.3 to obtain 
B 1 <c\T 1/8 \\Sf^\M\ 1 ^U\\\AQ x S\ 
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similarly 

^2 < ^(^) ll^ll ll^ll -K =L/8 ||^|| 3 / 4 ||^|| 

and 

2A\V8_ ii9i _ / 2 \V8, 



5 3 <c(^) \\Sf\AQ x 8\ + c(£) \\6\\^\AQx5f/\ 
It follows that 

3 

Y,B t = Ji + J 2 

i=i 

where by Theorem 2.3 with p = q = 2 we have 

'2A\V8 



* = *(^) 11*11(11*11 + 2||^||)|^| 

- t(^) 1/4||5||2(I1511 + 2]m)2 + ^ lAQx61 



2 



and Theorem 2.3 with p = 8/3 and q = 8/5 we have 

- 2 

a7 



^=cfA) 1/8 || 5 ||3/4 ( || 5 || +2 || t/ || ) | A g A5 |5/4 



< ^3/ 1 \V3 

~ 2 



(^) 13||5||2(l|5|l+2||t/|l)8/3 + i WA512 - 



Here, C 2 and C 3 are as defined in (3.10). It follows that 

±\\Q x 5\\ 2 + v\AQ x 5\ 2 < J3 + J4 

where by Theorem 3.8 

/ A \ 1/4 
J 3 = ^(^J ||5|| 2 (||5||+2||t/||) 2 

<^(^) 1/ V(yi|V^«)(||^ n )||e^)/ 2 + 2||t/||) ; 



z/ 4 Ai 
A \V4 



and 



<C 2 (-£-) \\5(t n )\\ 2 e^ t - t ^(R 1 / 2 e^ t - t ^/ 2 + 2K 1 / 2 y 



>/4 = C3(^-) 1/3 ||5|| 2 (||5||+2||t/||) 8 / 3 

1 \ 1/3 H^ \||2 „B(t-t n )(\\ su \\\Mt-tr,)/2 , oil ttII\ 8 / 3 



- C3 (^a7) ll^)l| 2 e^-*")(||5(t n )||e^- t ")/ 2 + 2||t/||) f 

< c 3 (4v-) I || V<«-«»> (W-*-)/ 2 + 2K 1 /2) 8 / 3 . 
V z/ Ai / 
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Then by Theorem 3.2 and the fact that P\8(t n ) = we have 

j t \\Qx5\\ 2 + v\\\Q>5\\ 2 < \\QxS(t n )fg(t - t n ) 

where g is the function defined in (3.9). Multiply by the integrating factor e ^ A (* _t ™) and 
integrate from t n to t to obtain 

\\Qx5(t)f <M(t-t n )\\Q x 5(t n )\\ 2 

where M is the function defined in (3.11). Let 7 = M(h). By our choice of t* and h we 
have 7 < 1. It follows that 

||5(t n+1 )|| 2 = \\Q x S(t n+1 )\\ 2 < lim \\Q x m\\ 2 

t/-t n+1 

< hm M(t-t n )||g A 5(t n )|| 2 = M(/i)||5(t n )|| 2 = 7||5(tn)|| 2 . 

Therefore ||5(t n+ i) || 2 < R, which completes the induction. 

To finish the proof note that under these hypothesis we have, in fact, proven 

\\5( tn )\\ 2 <rR- 

The proof now finishes as in (2.14). □ 

Corollary 3.10. Given any t* > there exists A large enough depending only on K, 
||5(to)||) v, ^ and t* such that for any h G (0, t*] the approximating solution u given by 
(1.1) with t n = hn converges to the reference solution U of (3.1) as t — > 00. 

Proof: First estimate g(r) from (3.9) as 



where 



Therefore 



M(r) < e-^ Ar (l + z/A x / 4 L 2 ^ e < W+" A )«d s ) 
7/3/3 + z/A 



where 

m(r) = (1 - L 2 A" 3 / 4 )e-^ + L 2 A- 3 / 4 e^/ 3 ) T . 
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Differentiating yields 

49(3 2 



m'(r) = -u\(l - L 2 A- 3 / 4 )e-^ T + ^A" 3 / 4 ^/ 3 ^ and 

3 



m"{r) = z/ 2 A 2 (l -L 2 A" 3 / 4 )e-^ + — — L 2 A 



3/4 p (7/3/3)r 

9 



Given t* > choose A large enough such that m'(t*) < 0. Clearly m" > for this 
value of A. Thus m! is a strictly increasing function. It follows that m'(s) < for s G [0, t*] 
and therefore 



M(h) < m(h) = 1 + / m'{s)ds < 1 



for /i G (0,t*]. Hence, the approximating solution u given by (1.1) with t n = hn converges 
to the reference solution U of (3.1) as t — > oo. □ 

Corollary 3.11. If we take r\ = and 



9 /3c|/M/3 



then the results of Theorem 3.9 hold where t* may be chosen independent of ||<5(to)||- 
Proof: When 77 = then \\S(t )\\ 2 = \\QU(t ) || 2 < A\ □ 

The bound in Corollary 3.11 would be comparable to the bound in Theorem 3.7 if the 
exponent of 8/3 were instead 2. The power 8/3 comes as a result of using the L 8 / 3 norm 
and the particular form of Agmon's inequality we have used in estimating the nonlinear 
term. Using the same proof technique with different interpolation inequalities in place of 
Theorem 3.4, this exponent could be reduced to as near 2 as one might like at the expense 
of increasing the constant c. 

Next we state without proof the analog of Corollary 2.7 for the incompressible two- 
dimensional Navier-Stokes equations. 

Corollary 3.12. Let t* be the bound given in Corollary 3.10. Suppose t n+ \ — t n < t* 
where t n — > 00 as n — > 00. Then the approximating solution u given by (1.1) converges to 
the reference solution U of (3.1) as t — > 00. 

We finish with the equivalent of Theorem 2.8 for the Lorenz system. This result is 
interesting because it shows that even if the approximating solution doesn't converge to 
the reference solution it is still bounded. This is striking because in the case of continuous 
data assimilation it is unknown whether the approximating solution is in general bounded 
or not. In particular, the comments before Theorem 3.5 in Olson and Titi [11] indicate 
that the approximate solution obtained by continuous data assimilation is not known to 
be bounded if it does not converge to the reference solution. 

Theorem 3.13. There exists M5 independent of h and depending only on \f\, O and v 
such that \\u(t)\\ 2 < M 5 /(l - e ~ uXlh ) for all t. 
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Proof: Multiply (3.1) by Au, integrate over O and apply the inequalities of Cauchy- 
Schwartz and Young to obtain 

It follows from Theorem 3.2 that 

at z/ 

The rest of the proof is similar to the proof of Theorem 2.8. □ 
Corollary 3.14. If 

. 9 f2cK 1 / 2 + c\\8(t )\\^ 
A > 



a; 



1/3 V 



7/ 



then there exists a bound depending only on K, ||5(io)||, v, O and A such that the 
approximating solution u obtained with t n = hn is bounded by M 6 for any h > 0. 

Proof: The proof is the same as the proof of Corollary 2.9 for the Lorenz equations. Note 
that if rj = as in Corollary 3.11 then A and M 6 may be chosen independent of ||5(to)||- d 



4. Concluding Remarks 

We have studied discrete data assimilation for the Lorenz system and the incompress- 
ible two-dimensional Navier-Stokes equations. Comparing the results we have obtained 
for discrete data assimilation to prior studies of continuous data assimilation we find the 
following. For the Lorenz system Pecora and Carroll [14] showed that continuous data 
assimilation of the X variable lead to convergence of the approximating solution to the 
reference solution as time tends to infinity. In Theorem 2.5 we provide a similar result 
for discrete data assimilation provided the update time interval h is sufficiently small. For 
the incompressible two-dimensional Navier-Stokes equations Olson and Titi [11] obtained 
conditions on the resolution parameter A under which the approximating solution con- 
verged to the reference solution. Theorem 3.9 states a similar condition on A that leads 
to convergence of the approximating solution to the reference solution provided h is suffi- 
ciently small. In Corollary 3.10 we also show that for any h > there is A large enough 
such that the approximating solution converges to the reference solution as time tends to 
infinity. Thus, discrete data assimilation has been shown to work under similar conditions 
as continuous data assimilation. 

A striking difference between discrete and continuous data assimilation is given by 
Theorem 2.8 and Theorem 3.13 which show that the approximating solution obtained by 
discrete data assimilation is bounded even if it doesn't converge to the reference solution. 
Although the approximating solution appears bounded in the case of continuous data 
assimilation for all numerical experiments performed to date, there does not yet exist an 
analytic proof of this property. Boundedness of an approximating solution that does not 
converge to the reference solution remains a conjecture for continuous data assimilation. 
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In this work we have found analytic bounds on the update time interval h for discrete 
data assimilation which guarantee that the approximating solution converges to the ref- 
erence solution. It is natural to compare our analytic bound with numerical simulation. 
For the Lorenz system Corollary 2.6 indicates that for t* ~ 0.000129 for the standard 
parameter values a = 10, b = 8/3 and r = 28. For these same parameter values Hayden [6] 
performed a numerical simulation of discrete data assimilation using the 150-digit-precision 
variable-step-size variable-order Taylor-method integrator [9]. In this work the maximum 
value for t* was found numerically to lie in the interval [0.175, 0.1875]. Thus, convergence 
of the approximating solution to the reference solution numerically occurs for values of h 
three orders of magnitude larger than those guaranteed by our analysis. 

For the incompressible two-dimensional Navier-Stokes equations bounds on A for con- 
tinuous data assimilation appear in [11] and [12]. That work shows that the approximating 
solution converges to the reference solution for values of A smaller than expected from the 
analysis. As similar techniques are used to treat the discrete assimilation in this paper, we 
expect our bounds on A to be similarly conservative and that the approximating solution 
obtained by discrete data assimilation will converge numerically for much smaller values 
of A than given by Theorem 3.9 and Corollary 3.11. We also expect our bounds on h to be 
conservative. A computational study of discrete data assimilation for the incompressible 
two-dimensional Navier-Stokes equations is currently in progress. 

We conclude by returning to the motivating problem of using satellite imaging data 
to initialize a weather forecasting model. In applications the value of h governing the 
time interval between consecutive observational measurements is generally much larger 
than the value of At used by the numerical integrator. Therefore, it is more realistic to 
treat the observational data as measurements occurring at a sequence of times t n as was 
done in this paper rather than as measurements occurring continuously in time. In the 
context of the Lorenz equations and the incompressible two-dimensional Navier-Stokes 
equations we have obtained similar theoretical results for discrete data assimilation as for 
continuous data assimilation. We hope that these results will shed light on the differences 
and similarities between discrete and continuous data assimilation and help guide future 
work in understanding more complicated problems. 
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